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. We introduce the spherical field formalism for free gauge fields. We discuss 

the structure of the spherical Hamiltonian for both general covariant gauge 
and radial gauge and point out several new features not present in the scalar 
field case. We then use the evolution equations to compute gauge-field and 
field-strength correlators. [PACS numbers: ll.10.Kk, 11.15.Tk, 11.15.-q] 
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Spherical field theory is a new non-perturbative method for studying quan- 
tum field theory. It was introduced in |j] to describe scalar boson interac- 
tions. The method utilizes the spherical partial wave expansion to reduce 
^ ! a general <i-dimensional Euclidean functional integral into a set of coupled 

one-dimensional, radial systems. The coupled one-dimensional systems are 
^| then converted to differential equations which then can be solved using stan- 

dard numerical methods.0 The extension of the spherical field method to 
fermionic systems was described in 0. In that analysis it was shown that 
the formalism avoids several difficulties which appear in the lattice treat- 
ment of fermions. These include fermion doubling, missing axial anomalies, 
and computational difficulties arising from internal fermion loops. This find- 
ing suggests that the spherical formalism could provide a useful method for 
studying gauge theories, especially those involving fermions. As a small but 
important initial step in this direction, we contribute the present work in 
which we introduce and discuss the spherical field method for free gauge 

1 Support provided by the National Science Foundation under Grant 5-22968 and the 
Deutsche Forschungsgemeinschaft 
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4 For free field theory these equations can be solved exactly, as we will demonstrate here 
for gauge fields. 
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fields. 

The basic formalism for spherical boson fields was described in 0. In this 
paper we will build on those results with most of our attention devoted to 
new features resulting from the intrinsic spin of the gauge field. We discuss 
the operator structure of the spherical Hamiltonian in detail, using two- 
dimensional Euclidean gauge fields as an explicit example. Like standard 
field theory gauge-fixing is essential in spherical field theory, and we have 
chosen to consider general covariant gauge and radial gaugef]. In each case 
we derive the spherical Hamiltonian and use the corresponding evolution 
equations to calculate the two-point correlators for the gauge field and the 
gauge-invariant field strength. Free gauge fields in higher dimensions can 
be described by a straightforward generalization of the methods presented 
here. The application of spherical field theory to non-perturbative interacting 
gauge systems and related issues are the subject of current research. 



2 Covariant gauge 



In this section we derive the spherical field Hamiltonian for general covariant 
gauge. We will use both polar and cartesian coordinates with the following 
conventions: 

t = (tcos6,tsm6) = (t\t 2 ) = (x,y). (1) 
In general covariant gauge the Euclidean functional integral is given by 



(UP A') exp 



dtL 



where 



L = J dOtl-lF^F^-^d^f] 



We can write the field strength F 12 as 

F 12 = h\(4- X + •£) (A* - iAy) - ( £ - ,JL) (A* + iA* 



>2i L e \dt + tde) A 



\/2-, 



-id (d_ 



A- 1 ] 



(2) 



(3) 



(4) 



'See H and references therein for a discussion of radial gauge. 
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where 

A x TiA y = ^2A ±l . (5) 
We now decompose A ±x into partial wavesQ 

A±l = 7T, E (6) 

n=0,±l,- 

Returning to our expression for the field strength, we have 

^ 12 = *A £ <?-(F?-F?), (7) 

n=0,±l,- 

where 

7-1+1 _ dA t-i n-1 a+1 ( q\ 

b n ~ —at r A n-l \P) 

F-l — 9A n+l I n+1 A-l (q\ 

We can also express the gauge-fixing term in terms of F^ 1 , 

d ^ = ^ E e M {^ + ^). (10) 

n=0,±l,- 

With these changes the Lagrangian, L, is 

I E ( F -n ~ F-n) - F- 1 ) - £ ]T + ^) + F- 1 ) . 

n=0,±l,- n=0,±l,- 

(11) 

In |l]] the spherical Hamiltonian for the scalar field was found by direct appli- 
cation of the Feynman-Kac formula. This is also possible here, but in view of 
the number of mixed terms (a result of the intrinsic spin degrees of freedom) 
we find it easier to use the method of canonical quantization. Let us define 
the conjugate momenta to the gauge fields, 



SL 



IT, 



+1 

n-1 a A +l 

5- 



I [(1 - h) F -n + ("I - h) F =n] (12) 



dt 

= = i [(-i - \)F? n + (i - . (is) 



U dt 



3 In our notation A^ 1 carries total spin quantum number nil. 
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Following through with the canonical quantization procedure, we find a 
Hamiltonian of the form 

n=0,±l,- 

where 

Hn — —Jt { n -n+l ~ "K-n-l) { n n-l ~ n n+l) (15) 
4f l"-n+l T l^B-l T Vlj f t /1 n-l 7I n-l ( /1 n+l''n+l- 

We obtain the corresponding Schrodinger time evolution generator by making 
the replacements 

A n l ^4\ ^-ajr. (16) 

We then find 

Hn = ~^ (sfc " " afe) (17) 



For the sake of future numerical calculations, it is convenient to re-express 
H in terms of real variables. Let us define 

Mr. 



X 



for n > and 



| ( 2 r7+l + 2 -n-l + Z n-1 + 2 -n+l) ( 18 ) 

\ { Z n+1 + 2 -n-l — Z n~l ~ 2 -rt+l) (19) 

^ { Z n+1 ~ Z -n-l ~ Z n-1 + Z -n+l) (20) 

H n = Yi ( Z n+1 ~~ Z -n-l + 2 n-l ~~ Z -n+l) (21) 

*° = ^ (^r 1 - *±i) ■ (23) 

The variables -u n and y n correspond with combinations of radially polarized 
gauge fields with total spin ±n, while v n and x n correspond with tangentially 
polarized gauge fields with total spin ±n. In terms of these new variables, 

H = H + J2 K, (24) 

n=l,2,- 
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where 

H o = -TtS* - ft& - \ ( u o£^ + x o£^) > ( 25 ) 



and 



H 'n ~ It (dvl + Ll) 2t (dul + dyl ) ( 26 ) 



ltd d d d 

Despite the somewhat complicated form, we know that the spectrum of H 
should resemble that of harmonic oscillators. To see this explicitly let us 
define the following ladder operators: 

7T+ L I 9 I d \ T j _ J_ / 2a-n+an d _i_ 2+n-an d \ „. _i_ „, \ 

(27) 

T/+ _J_ [ -2a-n+an d _i_ 2-n+an 9 , \ T/- _ _J_ / _9 8 \ 

V n V2\ 4(n-l) 9n„ ^ 4(n-l) 9*> n ^ "™ ^ ' V n ^ \y dUn dVn J , 

(28) 

y+ 1_ | 8 I 8 \ y— _1_ ( 2a—n+an d , 2+n-an d , , \ 

^" ~~ v 7 ! fe ny / ' ^« ~ ^ ^ 4(n+l) % n ~T 4(n+l) fe n ^ n Xn J ' 

(29) 

Y+ — _L | -2a-n+an 9 , 2-n+an g , _ \ y- _ J_ / _d d \ 

(30) 

for n > 1 and 

U o=~^ U o=^A + ^o, (31) 
^o + = -^a§7' ^ = ^afc + ^o. (32) 
We have constructed these operators so that 

[Un,U+] = [V-, V+] = [X~,X:] = [Y-,Y+] = 1 (33) 
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for n > 1 and 

[U -,U+] = [X ",X + ] =1. (34) 

All other commutators involving these operators vanish. We can now rewrite 
H and H' n (for n > 1) as 

H = \U+U, + iX+X " + f , (35) 



ttI _ n+1 TJ+TJ- _i_ i n+1 y+ V~ _|_ ItiV+V - _|_ 

n t n n ' t n n ' t n n ' t n n ' 



2n+2 
t ' 



(36) 



We now see that Hq is equivalent to two harmonic oscillators with level- 
spacing -, and H' n is equivalent to four harmonic oscillators, two with spacing 

and two with spacing It is important to note that H[, however, is 
quite different. The s-wave configurations are independent of 8, and so H[ 
does not have terms of the form z^ 1 -^! and zT. 1 ^t- Furthermore the gauge 

fields are massless and so H[ depends only on the derivatives of z^ 1 and 
Zq 1 . One consequence of this is that the spectrum of H[ is continuous, and 
calculations involving H[ and relevant boundary conditions must be done 
carefully. We will discuss these issues later in our analysis of correlation 
functions. 

In preparation for our calculation of gauge-field correlators, let us now 
couple an external source J to the gauge field, 



dJ9t 



i F v F v-i( d . A iy + A.j 



(37) 



The new spherical field Hamiltonian is then 

n=0,±l,- 

As noted in M the vacuum persistence amplitude is given by 



(38) 



Z(J) a lim (6|Texp 



dtH{J) 



\a) 



(39) 



where \a) and \b) are any states satisfying certain criteria. These criteria are 
that | a) is constant with respect to the s-wave variables Zq , z^ 1 and has 
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non-zero overlap with the ground state of H as t — > + , and \b) has non-zero 
overlap with the ground state of H as t — > oo.f\ Because the spectrum of 
H[ is continuous, in numerical computations it is useful to include a small 
regulating mass, /i, for the gauge fields and then take the limit 11 — > 0. 

3 Radial gauge 

We now derive the spherical field Hamiltonian for radial gauge. We take the 
gauge-fixing reference point, t , to be the origin, 



(f- t ) ■ A = t- A = 0. 



(40) 



We expect this gauge-fixing scheme to be convenient in spherical field theory 
calculations for several reasons. One is that non-abelian ghost fields in radial 
gauge decouple, as they do in axial gauge. In contrast with axial gauge, 
however, radial gauge also preserves rotational symmetry. As we will see, the 
spherical Hamiltonian and correlation functions in radial gauge are relatively 
simple. Since 



?.A=j*[e"A +1 + e-»A- 1 ], 
we can impose the gauge-fixing condition by setting 

A' 1 = -e 2W A +1 . 
With this constraint we express the field strength as 



f u = ^ J2 eine 



dA 
dt 



n-l I 1 A+l 
3* ~T t^n-l 



n=0,±l,- 

The radial-gauge Lagrangian is then 



L = J d6t[-\F 12 F 12 ] 



dt t-n-1 



9^1 | 14+1 
at t^-n-l 



n=0,±l, 



(41) 



(42) 



(43) 



(44) 



7 One caveat here is that |a) must lie in a function space over which the spectrum of H 
is bounded below. 
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We again follow the canonical quantization procedure. The conjugate mo- 
menta to the gauge fields are 



SL 



IT, 



+ 1 

n-1 BA+ 1 

at 



2t 



8A- 1 



dt 



-A 4. I A +1 
+ t A -n-\ 



and the radial-gauge Hamiltonian has the form 

ra=0,±l,- 

In the Schrodinger language the Hamiltonian becomes 



(45) 



(46) 



r — \^ \i a a a 

4tdz+\ dz +1 ' t Z n-l^~^~ 

' * L n— 1 — n— 1 n— 1 



n=0,±l, 

As before we now define real variables, 



n ~ y/2 V n-1 ~ ^-n-iy 
l Z n-l ^-n-lj 



for n > and 



(47) 



(48) 
(49) 



Xq = i ■ z^\. 



(50) 



Our Hamiltonian can be re-expressed as 



H = H + ^2 H 'w 

71=1,2,- 



(51) 



where 



ff - i a 2 i„, a 



(52) 



it 



a 2 , a 2 



r — + ?; — 



(53) 



Let us now define ladder operators 

X+ — 1JL v- 

n 2 dx„ ' n 



AJt Or 

2 cte„ 



(54) 



8 



for n > and 



™ 2du n ' 2di>„ 



- 2f r. 



(55) 



for n > 0. These ladder operators satisfy the relations 

[X~,X+] = [V-,V+] =1, (56) 
while all other commutators vanish. In terms of these operators we have 



(57) 



\XnX n 



iv + v~ + 



(58) 



We note the radial gauge constraint has removed the continuous spectrum 
from the s-wave sector, and the spectrum of H is purely discrete. Further- 
more the splitting between energy levels of H' n is independent of ri.| 
We now couple an external source J to the gauge field, 



d9t 



12 T712 



2 L 



+ A-J 



The new Hamiltonian is then 



(59) 



(60) 



n=0,±l,- 

The vacuum persistence amplitude in radial gauge is given by 



Z(J) a lim (6|Texp 

/ ■ >Q+ 

'rnin T u 



dtH(J) 



^min 



(61) 



where \a) has non-zero overlap with the ground state of H as t — > + , and 
I b) has non-zero overlap with the ground state of H as t — > 00. Since the 
level-spacing of H[ diverges as t — > + only the ground state projection of 
H[ at £ = contributes, and |a) no longer needs to be constant with respect 
to the s-wave variables z^ 1 . This a rather important point since, as we recall 
from m, the constant s-wave boundary condition at t — follows from the 
fact that the value of field at the origin is not constrained. This is however 
not true in radial gauge as a result of the gauge- fixing constraint. 

8 The reason is that in radial gauge the gauge fields are tangentially polarized, and so 
purely tangential excitations in two dimensions do not contribute to the field strength F 12 . 
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4 Gauge-field correlators 

In this section we calculate two-point gauge-field correlators using the spher- 
ical field formalism. This calculation can be done in several ways, including 
numerically. For future applications, however, it is useful to have exact 
expressions for use in perturbative calculations (e.g., for evaluating countert- 
erms) . Here we will obtain results by decomposing the fields as a combination 
of the ladder operators. Let us start with radial gauge. We have 



<-l = iT2 ( tX n + iX n ~ Vn ~ V~) (62) 



Z+ -n-l = ^2 ( iX n + iX n + + V~) (63) 



for n > and 



z. 



+i 



i 1 = | W + *o) ■ (64) 



We would like to calculate the correlation function, 

/r-Vi(M') = (01^(0^(0 \0) rad . (65) 
By angular momentum conservation f™\ n ,_ x vanishes unless n' = —n. Also 

/iViM^rVif^). (66) 

and so without loss of generality it suffices to consider n-1 for n > 0. 

For typographical convenience let us define 

(iF\ iG\ ) = 6(t- t') ( b \ U (°°' t ) FU ( t ' t ') GU ( t '> )\ a ) _j_ Q( t > _ A (b\U(oo,t')GU(t>,t)FU(tfi)\a) 
\\ ft 1 ft' I v / (6|f7(oo,0)|a) V > (6|f7(oo,0)|a) 

(67) 



where 



Ufah) =Texp 



*2 



dtH 



(68) 



For n 7^ 0, we have 



/r-l-n-i = I <{^n + + iX~ - V+ - V~} t {iX+ + iX~ + V+ + V~} t ) . 

(69) 
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Using the commutation properties of these ladder operators with H(t), we 
find 



For the special case n = 0, we obtain 



(70) 



f-iU = I <{^o + + iXo} t {^o + + iXo} v ) = ~\ W - t)7 + W - 0?] • 

(71) 

These correlation functions are in agreement with results we obtain by de- 
composing the following expression |§ into partial waves: 



(01^(2)^(0)10) 



racl 



t- lim 



•Itt 



£^0 + 



(72) 

(*-yf _ r 1 j a n,5 1^ (gg-g)!±g _ r 1 w+ „< (^-^) 2 + £ 



^ l og 1£Z|L _ df £ ds X J log i££^+£ _ 9 | £ dt y i log _^ 

+«9f £ rf S jj f • y log 



The length scale L is used to render the argument of the logarithm dimension- 
less. Its purpose, however, is only cosmetic since the gauge-field correlator is 
not infrared divergent in radial gauge and the dependence on L cancels.^ 

These same methods can be applied to gauge-field correlators in covariant 
gauge. We find 

<0| A+\{t)AtU{t') |0)_ = <0| AZUMA-Ut') |0)_ (73) 

= 2^,0 P(*-0t + o(t'-tH] 



s=i[0(i-f)Vi + fl(f-%,i]- 



and for n^l, 

(01^,(^)^(010) 



q+l 
cov 4|n— 1| 



9{t-t>) (i) ln - ll + 9(t'-t) (1) 



|n-l| 



In covariant gauge the s-wave correlator is infrared divergent, 
(0| A+\t)A^(t>) \0} cov = -s±l U - log | + 0(f - 1) log £ 



(74) 



(75) 



9 The radial gauge constraint A 1 = —e 2l9 A +1 pairs s-wave configurations with d-wave 
configurations. Since the d-wave is not infrared divergent, neither is the s-wave. 
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where logL is infinite. This divergence is specific to two-dimensional gauge 
fields and does not occur in higher dimensions. If we include a regulating 
gauge-field mass, /i, we find that L scales as 1/fi as fi — > 0. These correlation 
functions are in agreement with the results we obtain by decomposing the 
following known expression into partial waves: 

(0\A i (a!)Al(if)\0) eou = -± 

5 Gauge- invariant correlators 

Let us now consider the two-point correlator of the gauge-invariant field 
strength F 12 . We can calculate the F 12 correlator by differentiating the 
gauge-field correlators calculated in the previous section, but it is instructive 
to redo the calculation by coupling a source to F 12 . This time we describe 
the calculation in detail for covariant gauge. The same calculation can be 
done for radial gauge using similar methods. Let us start by quoting the 
result we expect. From free field theory we know 

(0| F 12 (t)F 12 (?) |0> = 5 2 (t- f). (77) 

The F 12 correlator has a simple local structure, a consequence of the fact 
that in two dimensions gauge fields can be decomposed into scalar and lon- 
gitudinal polarizations (borrowing Minkowski space terminology). There are 
no transverse polarizations to produce non-local contributions to the F 12 
correlator. We will see this happen explicitly in the calculations to follow. 
The two-point correlator for partial waves of F 12 is given by 

(0| F^ 2 (t)F r ]?(t') |0) = ^ J d6d6'e- ind e- in ' ' (0| F 12 (t)F 12 (f ) |0> , (78) 

and we deduce 

<0| Fl\t)Fl 2 {t') |0> = |^(t - if). (79) 

Let us now reproduce this result using the spherical field method. We return 
to the covariant gauge Lagrangian and couple a source /C to F 12 , 

L — f d6t [-\F 12 F 12 - ±(d t A l ? + F 12 K] . (80) 
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The conjugate momenta are now 



7T 



-1 



t 



7T. 



re— 1 2 

-1 _ t 
re+1 2 



(1 - + (-1 - i)Fll - iv^C-n" 

(-1 - + (1 - + iV2/C-» 



and the new Hamiltonian is 

H{K) 

where 



n=0,±l,- 



H n {K) = H n (0) + -^/C n 7T n _ 1 - -^/C n 7T n+1 - ^K- n K n . 

In the Schrodinger language we have 

H n (K) = H n (0) + j=K n {-^^ + i^) - |/C_ n /C n 
for n^O and 



H ()C) = H Q (0) - )C £- o - IK )C . 
The vacuum persistence amplitude is 



Z(/C) = (6|Texp 



dtH{K) 



| a) 



(81) 
(82) 

(83) 

(84) 
(85) 
(86) 

(87) 



and we evaluate the F^ 2 correlator by functional differentiation with respect 
to /C, 



(o| i^(f)i#(0 |o> 
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Z{0) tdK- n (t) t'dK_ n ,{t>) 



K=0 



It is clear from angular momentum conservation that this correlator vanishes 
unless n = —n', and so it suffices to compute 



(0\F^(t)F^ n (t') |0>. 



(89) 



Since (|89|) is symmetric under the interchange n, t «-> —n,t', we can also 
restrict n > 0. Differentiating with respect to the sources, we obtain, for 
n > 0, 



(o|^(t)^(t')|o) 



i 



— + i— 



dv 



(90) 
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When n > 1 we can write + ijf- as a linear combination of £/+ K , 
X+, and Y~. These ladder operators are, however, acting in four different 
spaces. The matrix element of the operator 



u(oc, t) (£- + i£.)u(t, t')(-£ + i£-)u(t f , o) 



(91) 



' dx n > \ ' A dv n dx rt 

from the ground state at t — to the ground state at t = oo vanishes.^ 
Consequently for n > 1 only the delta function contributes to the correlation 
function in (p0|). The same arguments apply for the case n = 0, and only the 
delta function contributes here as well. 

We now turn to the special case n = 1. The relevant part of the Hamil- 
tonian is 



\_ I a 2 , a 2 

2t \ dv^ dxj 



a 2 , a 2 

2t V du 2 dy'f 



(92) 



The combinations U\ 
of the s-wave variables Zq and z, 



v\ and X\ — y\ correspond with linear combinations 
1 . The initial configuration |a) at t = is 



constant with respect to z^ 1 and z$ , and therefore constant with respect to 
ui — v\ and x 1 — y\. We note that when H[ and/or ±^j- + acts upon 
|a), the result is again a state constant in ui — v± and x% — y\. It therefore 
suffices to compute the correlator restricted to the subspace which is constant 
in U\ — v i and X\ — y\. In this space H[ has the form 



Q + l 

8* 



a I a 

dui dvi 



q+l 
8t 



9 I a 

dxi dyi 



(93) 



dyi 



Comparing with (p5|), we see that this is analogous with the previous case 
for n = 0. We again find the result 



<0|F^(t)F«(O|0) = itf(t-O 



(94) 



6 Summary 

In this work we applied the methods of spherical field theory to free gauge 
fields. We analyzed two dimensional gauge fields in general covariant gauge 

10 As noted before, this is due to the fact that in two dimensions there are no transverse 
polarizations. 
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and radial gauge. In the process we have discussed several new aspects which 
resulted from the spin degrees of freedom as well as the masslessness of the 
gauge field. As we have seen, polarization mixing complicates the structure 
of the spherical field Hamiltonian. Nevertheless in radial gauge we were able 
to decompose the spherical field Hamiltonian as a sum of harmonic oscilla- 
tors. We did the same for covariant gauge, but found that the s-wave part of 
the Hamiltonian has continuous spectrum. In relation to these differences, 
we also discussed issues regarding the s-wave boundary condition at t — 0. 
We then used the spherical field evolution equations to calculate two-point 
correlators for the gauge fields and field-strength tensors F 12 . Our presen- 
tation here is intended as a first introduction to the application of spherical 
field methods to gauge theories. Free gauge fields in higher dimensions can 
be treated by a straightforward generalization of these methods. Interacting 
gauge systems, however, include many interesting theoretical and computa- 
tional issues not discussed here, and these are the subject of active research. 
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